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We derive an asymptotic expansion for a Wiener-Hopf determinant arising in the problem of 
counting one-dimensional free fermions on a line segment at zero temperature. This expansion 
is an extension of the result in the theory of Toeplitz and Wiener-Hopf determinants known as 
the generalized Fisher-Hartwig conjecture. The coefficients of this expansion are conjectured to 
obey certain periodicity relations, which renders the expansion explicitly periodic in the "counting 
parameter". We present two methods to calculate these coefficients and verify the periodicity 
relations order by order: the matrix Riemann-Hilbert problem and the Painleve V equation. We 
show that the expansion coefficients are polynomials in the counting parameter and list explicitly 
first several coefficients. 
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I. INTRODUCTION AND MOTIVATION 

Toeplitz determinants are determinants of matrices whose elements depend only on the difference of the matrix 
indices: 



They occur in many topics of theoretical physics: statistical physics [3 [2], random-matrix theory [3i 4 , full counting 
statistics of fermionic systems [5 , non-equilibrium bosonization [5], etc. Typically, in physical applications, one is 
interested in the behaviour of a Toeplitz determinant ([T]) as the matrix size N tends to infinity. For rapidly decaying 
matrix elements a^^j, the leading exponential dependence oi Dn on N can be easily understood on physical grounds: 
the coefficient in the exponent is given by the average logarithm of the "symbol" (the Fourier transform of ai_j) of 
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This resuh is known as the strong Szego theorem ^i8j [the theorem also gives the coefficient A in terms of (j{k)]. li the 
matrix elements ai-j decay slowly (or, equivalently, if the symbol cr[k) has singularities), the exponential dependence 
([2]) is complemented by power-law prefactors: the result known as the Fisher-Hartwig formula [71I5HT^. 

A particularly interesting feature of the Fisher-Hartwig formula comes from the ambiguity in choosing the branch 
of the logarithm in Eq. ([2| in the case of a singular a{k). As a result, one obtains multiple branches of the asymptotic 
dependence of Z^at on TV, and the traditional Fisher-Hartwig formula prescribes selecting the leading one (having 
the prefactor with the largest power of A^). If the Toeplitz matrix depends on a parameter, the leading branch may 
switch as a function of this parameter: in this case, the asymptotic behavior of Dn depends on it nonanlytically. At 
the switching point, two branches are equally relevant, and the correct asymptotic behavior is given in this case by 
a simple sum of these two branches. This prescription known as the generalized Fisher-Hartwig conjecture was 
recently proven in Ref. [T?l 

In recent literature, it was conjectured that subleading branches of the Fisher-Hartwig formula do not need to 
be discarded, but they provide an accurate description of subleading terms in an asymptotic expansion of Dn as N 
tends to infinity |15H18) . Moreover, each of the Fisher-Hartwig branches ([2]) may, in turn, be improved by including 
corrections as a usual power series in 1/N. It was further conjectured that a full asymptotic series for Dpf may be 
obtained as the sum of all Fisher-Hartwig branches, in which all terms in the 1/A^ expansion are kept [T^ [50]. A 
particular case of this conjecture was also proposed in Ref. |5] for the Toeplitz determinant describing the full counting 
statistics of one-dimensional free fermions. It was verified numerically that, in this example, the first subleading 1/N 
terms in the leading and subleading Fisher-Hartwig branches reproduce several terms in an asymptotic expansion of 
Dn. 

In this work, we support this conjecture by an explicit calculation for the problem of free fermions on a continuous 
line, which is a limiting case of the lattice model considered in Ref. [51 In this limit, N is replaced by a continuous 
parameter x, the Toeplitz determinant becomes a Fredholm determinant (more specifically, a Wiener-Hopf determi- 
nant with a piecewise constant symbol and we may use methods developed for the latter [2 HH [Ml US] • We 
derive a complete asymptotic expansion consistent with the above-mentioned genrealization of the Fisher-Hartwig 
conjecture. A proof of this new conjecture is still missing: it amounts to certain "periodicity relations" on the coeffi- 
cients of this expansion. However, we present a systematic algorithm for calculating the coefficients to an arbitrarily 
high order in 1/x, which allows us to verify these periodicity relations order by order. We have verified the periodicity 
relations up to the 15th order in and conjecture that they hold to all orders. 



II. FORMULATION OF THE PROBLEM 



We begin our definitions with a description of a relevant physical problem of full counting statistics of free fermions 
on a segment of an infinite line. Consider free fermions in one dimension at zero temperature. Their multi-particle 
state is characterized by a single parameter: the wave vector kp such that all states with wave vectors |fc| < kp are 
filled and all states with > kp are empty. We are interested in the expectation value (exp[27rzK(5]), where Q is the 
operator of the number of particles on a given line segment of length L and k is an auxiliary "counting parameter" 
[51 |23H15] . This expectation value may be re-expressed as a determinant of a single-particle operator |51 [5SH2H] : 

(exp[27riKQ]) = det [l -I- np (e^'^^^O - l)] = dot [l + npQ (e^"" - l)] . (3) 

Here Q in the right-hand side is the single-particle projector on the line segment in real space, 

^ fl ifO<,<L, 
I otherwise 

(q is the coordinate on the line) and np is the projector on the occupied states in the Fourier space, 

Jl if|fc|<fcj., 

np — < (5) 
otherwise . 
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The determinant (|3| is of the Wiener-Hopf type pTl 122) . It may be understood as a Fredholm determinant, 

det[l + npQ{e'^''' ^l)]=y^^—— dq, . . . dqi det K(g,-<z,)], (6) 

Jo Jq 1<»J<' 



where np{qi — qj) is the Fomier transform of np defined by Eq. ([s]): 

sin(fcF(7) 

np{q) = . (7) 

irq 

Obviously, the determinant ^ depends on kp and L only via their product kpL, and therefore we may define 

x{k,x) =det[l + npQ{e'^'''^ -I)] , x^kpL. (8) 

Alternatively, we may define the same function by discretizing the space coordinate q and then taking the 

continuous limit (as in Ref. [Sj. Namely, we may first consider the Toeplitz determinant on a lattice: 

DN{K,kp)^ det [5,j + np{i - ]) {e"^'^ - I)] (9) 

l<i,j< 

[here np(i — j) is defined by the same expression ([7|] and then define the function x('*,a^) as the limit 

x{k,x)= lim Dn{k,x/N) (10) 

This type of definition is refered to as a "double-scaling limit" in Ref. TH 

We have thus defined a function of two variables xl'^j ^)- At a given x, it is periodic in k with period one. From the 
expansion ^ it follows that xi^^j ^) is an entire function of k at any fixed value of x. The main result of this paper is 
a conjecture of an exphcitly periodic asymptotic expansion for x(K,a;) at a; — )■ oo, Eqs. (11) and (12). As announced 
in the introduction, the leading asymptotics of xi'^i a;) at a; — )■ oo is discontinuous in k (at points ^{k) — 1 + 1/2 with 
integer /), and our asymptotic expansion describes in full detail the development of this discontinuity related to the 
switching of Fisher-Hartwig branches. 



III. MAIN RESULTS 



In Section IV below, we derive the asymptotic expansion (34) for the function x('*i x)- This form of the asymptotic 
expansion is proven, provided ^{k) 7^ Z-l-1/2 for any integer and the coefficients Lj^n, x) are computable, as Laurent 
series in l/x, iteratively order by order. 



Furthermore, we conjecture that the coefficients Lj{K,x) obey the "periodicity relations" (35 1, which brings the 
expansion (34) to an explicitly periodic form. Under this assumption (which we verified to many orders in l/x), the 
expansion ( 34 1 may be brought to the form 



+00 



x{k,x) = ^ x*{ii+j,x) , 



(11) 



where 



X*iK,x) = exp 



2iKX — 2k^ \nx + C{k) + /ji('t) (ix) 



n=l 



(12) 



This asymptotic expansion agrees with the general conjecture proposed in Refs. [13 [20! and with a more explicit 
formula conjectured in Ref. i5l 

The sum in Eq. ( 11 ) corresponds to adding together all different Fisher-Hartwig branches, and the sum in Eq. ( 12 ) 
includes all (1/a;)" corrections within a given branch. The coefficient C{k) is given by [51 155] : 

C(k) = 21n|G'(l + k)G(1 - k)\ - 2k^Iii2, (13) 
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where G{z) is the Barnes G function [55]. The coefficients are polynomials in k with real rational coefficients, 

and they are odd/even in k at odd/even n, respectively. Moreover, the lowest power of k in /„(k) is 3 or 4 (for n odd 
or even, respectively) [30] • The first several coefficients are [31) : 



/i(« 

/3(« 
/4(k 
/5(k 
/6(k 
/7(aJ 



2" 
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K 
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(14) 



12k^ + 
1931 
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24 
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16' 



8263 7 2155 . 45 , 

-K + — rr— K + — K . 



224 16 32 56 

In Sections IV and [v| we give algorithms for calculating the coefficients fn{n) order by order up to an arbitrary large 



In this work, we prove neither the expansion (11)-(12| nor even its weaker form (34) at points 5ft(K) = I + 1/2. 
However we conjecture that it is also valid there, so that the expansion (11)-(12) is in fact a uniform asymptotic 
expansion on any compact subset of k. 



IV. DERIVATION USING THE MATRIX RIEMANN-HILBERT PROBLEM 



In this section we show how the expansion (111, (12 1 can be derived by means of the asymptotic solution of the 



associated Riemann-Hilbert problem. The relationship between the the matrix Riemann-Hilbert problem, integrable 
partial differential equations and Fredholm determinants of integrable kernels, in particular the sine kernel, was 
exploited in different contexts such as classical inverse scattering problem, random matrix theory and quantum 
integrable systems [MJ 1321 - 135] . The asymptotic solution of the matrix Riemann-Hilbert problem for the generalized 
sine kernel with Hartwig-Fisher singularities was developed in Refs.|T31[Tni[3ni In this section we follow the notations 
of Ref. EU In that work, the logarithmic derivative of the function %(«;, x) was expressed in terms of a solution of a 
certain Riemann-Hilbert matrix problem, and then this problem was solved to the leading order in x. In this section, 
we solve the same Riemann-Hilbert matrix problem in terms of a series in 1 /x whose coefficients may be calculated 
iteratively, order by order. 

We do not repeat here the full argument of Ref. (391 but start with formulating their result relevant for our calculation. 
We refer the reader to the original paper for its derivation. The logarithmic derivative of x(k, x) was expressed there 



dx 



lnx('*, x) — 2iK 



i lim k [5oo(fc)ii - 1] , 



(15) 



where >S'oo(fc)ii is the upper-left-corner matrix element of the 2x2 matrix Soo{k), which solves the Riemann-Hilbert 
problem (in the complex variable k) formulated below. 

Define regions T>r and as Refc > /c^ and Refc < k^, respectively, where the region boundaries are chosen as 
— l<fcL<fcij<l (see Fig. [T i). This choice of regions differs slightly from Ref. '35', where the regions were chosen 
as discs around k = ±1: this difference does not change anything in our calculation, but simplifies the discussion of 
branch cuts. In these regions, we define the two matrix-valued functions: 



bR[x{k + l)]-"- 



V x{k-l) 



—ix{k 

mi 



1)] 

K, ~ix{k — 1)] 



aR[x{k + 1)]' 
xik~ 1) 



-2k 



-*[1 + K,ix{k - 1)] 



^[K,ix{k - 1)] 



(16) 
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FIG. 1: (a) Initial choice of domains for the matrix Riemann-Hilbert problem (21 1. The dashed lines denote th e b ranch cuts 



of 6lfl(fc) in T>R and of O^ik) in Vl- The arrows at the boundaries 92?^,^ denote their orientation in the integral (301. (b) The 
choice of domains used in Ref. [39 and in our formula (331. 



in the region and 



-Xk) 



^[k, ixik + 1)] 



aL[x{k-l)Y\~ 



x{k + l) 

*[1 + K, -ia;(fc + 1)] ^[-n,-ix{k + l)\ 



\ x{k + l) 

in the region T)]^. Here the constants and bji are the same as in Ref. "39, 



^l- K,ix{k+i)] ^ 
J 



(17) 



an = - 

the coefficients and &i are defined as 

bL - 



r^{-K) ' r2(K)sin2(7rK) 

'^^7 r, ai = -bflsm (ttk) , 

sm (ttk) 



and we have introduced a shorthand notation 



(18) 



(19) 



(20) 



where ^'(a, l;w) is the Tricomi fmiction j40, • The functions ^(a, l;?!") and "^(a^w) are defined to have branch-cut 
discontinuities along the negative real axis w G M_. Then the matrices OR{k) and 0i,(A;) have discontinuities along 
the Hues Refc = ±1, respectively. Note however that those discontinuities decay exponentially in \x{k =F 1)1, so that 
9R{k) and B^ik) have regular asymptotic expansions in inverse powers of x{k =F 1) (see Ref. l39l for more detail). We 
now formulate the Riemann-Hilbert problem for a matrix S{k): 



(i) S{k) is analytic in C\{dT> r\J dV r) 

(ii) S'oo(fc) = SR,L{k)eRx at A: e dVR^L ] 

(iii) Soo{k) I as fc — > cx) , 



(21) 



where we denote by Sr^l,oo the restrictions of S{k) onto Vr, 2?^, and C\{T>rUT>l), respectively. The matrix Soo{k) 
obtained as a solution to this Riemann-Hilbert problem can be then used to find the logarithmic derivative of x(k, x), 
according to Eq. (15). This relation is exact and is a minor reformulation of results obtained in Sections 3 and 4 of 
Ref. [31 

We can now calculate the asymptotic expansion for Soo{k) by expanding Or^l in powers of l/x and matching the 
matrix at the contours dT>R L order by order. In this expansion, we will treat all powers of and of e" as terms 
of order zero with respect to x: in other words, we will collect together all terms with the same integer powers of 
x, while letting the coefficients to depend on x'^ and e^^ . As we shall see below, such a method indeed produces an 
asymptotic expansion for Sco{k) within the interval |5i(K)| < 1/2. 

We start with the (asymptotic) expansion 



§(a, w) 



^ Pn{a) 

n=0 



where 



Pn{a) 



TJ,' 



(a + n- 1)2 



(22) 



6 



to expand 



1«« +1 



i(2) 



where 



and 



o(n) 



(fc - 1)» V *"-^&flP„-i(l - [x{k + 1)]"" 



(-z)"-iaflP„_i(l + K) [a;(A: + l)]" 



(fc + 1)« V «""'&LPn-i(l + - 1)] 



(-z)"-iaiP„_i(l-K) [x(fc-l)]' 



(23) 



(24) 



(25) 



(note that fl^^j^ depend on a; themselves, but only "weakly", via x^^'*). 

Now the Riemann-Hilbert problem (21 ) may be solved iteratively (order by order) in terms of an expansion 



5(fc) - / + -^W (fc) + ^5(2)(fc) + 



(26) 



where S'*^"-'(fc) are polynomials in x^^". Note that this method slightly differs from the approach used in Ref. [39l 
where the ansatz (3.53) allowed to partly resum the series (26). We denote the functions 5*^"' in the three domains 



Vr, Vl, and C\(I?fl U Vl) by 5*^^^ s'^l^^, and 5*^^ respectively. To the first order, we find 



while the general equation at the n-th order is 

StHk) = ^^^(fc) + + ■ • • + 4'L(fc)<Z'^(fc) + at k G dVn,L . 



(27) 



(28) 



Using the analyticity of S'"j^\ and S^' in the domains 1?^, V^, and C\(I?fl UI?l), respectively, and the boundary 



An) 



solving Eq. (27), we find 



condition S^\k) at A; — > oo, we can solve these equations by the Cauchy integral formula. At the first order, 

dk' 1 ^ 

Q-Pft ' ~ JdVL ^■""^ k' — k 



S^^\k) 



dk' 1 
27rj k' — k 



ik' 



?W(A:') 



(29) 



and, more generally, at the n-th order, the solution to Eq. (28) reads 

/ dk' 1 
Ig-f)^ 2-111 k' — k 



a=R,L 



s^^Hk')ol^''\k') + e^^\k') 



(30) 



These formulas produce the components S^\k), S^j^\k), and S^j^\k) depending on the location of the point k. At 
this stage of the calculation, it is technically convenient to deform the integration contours dV^ into the boudaries 
dVa of some nonoverlapping discs Vn and centered at A: = 1 and fc = — 1, respectively, as in Ref.|Ml(See Fig.[l]D). 
This deformation is allowed, since the ingtegrations converge rapidly at infinity and the cuts of the matrix functions 
(^R.hik) disappear from the asymptotic expansion. 



The formula (30), in principle, solves our problem: knowing the expansions (24) an d (25), we iteratively calculate 
S^'^\k) from (30) and then extract the logarithmic derivative of x(k,x) using (15). As a result, we obtain the 
asymptotic series [41] 



d_ 
dx 



lnx(K, x) — 2iK + In{K, ; 



rs=l 



where 



/„(k,x) — ^ ^ Rn,m (^) 



(31) 



(32) 
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is the contribution from S^ai^k). The form (32) of the term /„(k,x) follows from examining the explicit expression 

dkg 1 
Ido' ' '' 



{Uj} {Uj} 

nj>l aj—B.,L 
J2j = i "3=" 



dki 



27rz ki — /s2 



1 



k 



s-l 



9t'\ki)...e(-;\k,) 



(33) 



which is obtained by an iterative application of Eq. ( 30 ) . Here the sum is taken over all integer partitions of n into the 



sum of rij, and, for each such partition, over choices of the left and right integration contour for each j . Furthermore, 

the integration contours are ordered in such a way that the contour dD^ lies inside dDa ^, if j > j'. Every integral 
in Eq. (33) can be easily calculated by residues. The coefficients i?„,m(K) can be easily extracted from this expression 



by selecting terms with a particular power of x. The following properties of these coefficients can be proven: 



1. The coefficients Rn.min) vanish for |r7i| > [n/2]. This means that the sum over m in Eq. (32) extends only from 
-[n/2] to 



(k). This follows from the symmetry 



V2]. 

2. The coefficients Rn.m{i^) have a definite parity: Rn,m{~ii) = (— l)"~^-Rn,-7n 
k ^ —k (see a detailed discussion of this symmetry in Section 3 of Ref . [5^ . 

3. The coefficients Rn,o{n) are polynomials in k with real rational coefficients multiplied by (— The smallest 
possible degree of k contained in this polynomial is 3 or 4, depending on the parity of n. This can be easily seen 
for every term in the sum (33) by using the identity anhji = ULbL = —k?. 



4. The first coefficient in the expansion (31) is Ri q{k) = —2k} [calculated directly using the formula (33)]. 



Note that the powers of x in the formal series (31)-(32) decay if and only if |5R(k)| < 1/2. Therefore our calculation 
produces an asymptotic expansion for [d /dx) lnx(K, x) only within this interval of values of k. 



To obtain the expansion of the function x(K,a;), we integrate and then exponentiate the expansion (31), (32) as a 
formal series in x~^ . In the resulting series, we collect together terms with the same oscillatory prefactor e^*^^. The 
result may be further written as 



+00 

E 

j = -oo 



exp [2i(K + j)x — 2{k + j)^ Inx] Lj{K, x) 



(34) 



where Lj (k, x) are Laurent series in x with coefficients depending on k. These coefficients may be expressed in terms of 
Rn,m and vice versa, modulo an overall numerical prefactor in all Lj(K,x), which is left undetermined [it corresponds 
to the integration constant of Eq. (31)]. 

It seems very plausible [and it was conjectured both in the Toeplitz (chain) and Wiener-Hopf (continuous) cases 
[51 1191 [^ ] that the expansion of the form (34) is explicitly periodic in k, namely 



Lj{K,x) = Lo{k + j,x) 



(35) 



We do not have a proof of this conjecture at the moment, but we have verified it analytically up to the order x in 



Eq. (34) using the technique based on the Painllcve V equation, see Section |v] We conjecture that the relations (35) 
hold to all orders. 



Finally, we observe that Lo{k,x) does not contain negative powers of x [and, under the periodicity conjecture (35), 
neither do any of Lj(K, x)] and, therefore may be formally written as 



Lo{k, x) — exp 



c^('*) + E/"(«)(*^)"" 



n=l 



(36) 



(we included the factors i " in the expansion to make the coefficients real). Under the "periodicity conjecture" (|35|, 
this immediately leads to the perioduc form of the expansion ( 11 )-( 12 ). 



The coefficient C (k) cannot be calculated by the method described, but it is known from other approaches [51 

and is given by Eq. (13). To calculate the coefficients /«('*), we relate them to Rn.m{n). In fact, it is sufficient 
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to consider only the coefficients i?„^o(K). By comparing the expansion (31), (32) wtih (111 - (131, one finds: 



where we denote 



and 









i?3,o(K) 






■R4,o(k) 








= -4»-4 [/4(«;) 




Re, o(k) 


= -5z-' [/5(k) 




-R7,o(k) 


= -6r^ 






-R8,o(k) 


= -ll-' 


/7(«) 






- exp [C(k 


+ 1)^ 





(37) 



A2 

A3 



A? 
2 



A2A1 



2C(k) 



A? 
6 



16 



A„>i = /„(k + 1) + /„(k - 1) - 2/„(k) . 



(38) 



(39) 



Note that i?„^o(K) with n > 4 contain cross terms (containing A„) arising from Fisher-Hartwig branches with j 7^ 
in Eq. ( 11 ). By expUcitly calculating Rnfiin) from Eq. (|33|) [we used a computer program to perform this calculation] 
and solving Eqs. (37), we arrive at the results ( [l4| . Wealso remark that, for large n, the use of the formula (33) is 
not practical for explicit calculations (the number of terms grows very rapidly with n), and Eq. (30) seems to be more 
efficient. 



We can now prove the properties of the coefficients /n(K) declared in Section III 



• The reality of /«(«) follows from the property 3 of Rnfl{K): indeed, if we choose (ix) as the expansion variable, 
then all the coefficients of the expansions become real. 

• The parity of the coefficients /„(k) follows from the parity of i?„.o('*) (property 2). One can also formulate this 
property as an invariance of the whole asymptotic series with respect to the simultaneous formal sign change of 



X and K in their integer powers, while transforming fractional powers of x as x 1—^ x 



-4k 



• One can also easily verify that the cross terms in relations (37), for any Rn,Oj ^-re always polynomials divisible 
by K^. Therefore the property that i?„^o is always divisible by also holds for /„(k). 



V. CALCULATION USING THE PAINLEVE V EQUATION 



An alternative way to obtain the asymptotic expansion (11), (12) is the use of the Painleve V equation. It was 
discovered in the seminal paper |43j (with a simpler version of the derivation presented later in Ref. that the 
Fredholm determinant ([s]) considered as a function of x satisfies an ordinary differential equation: the Painleve V 
equation in the Jimbo-Miwa form, 



ixa"f + 4{xa' - a) {xa' -a + (a')') = . 
Here prime means the derivative with respect to x and 

d 

a{K,x) = X— lnxiK,x) . 



(40) 



(41) 



Remarkably the parameter k does not enter the equation ( 40 1 itself but defines its solution through the boundary 
condition: 



cr(K, x) 



1 



1 



x^ + Oix^) 



as X ^ 



(42) 
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[the same expansion can also be obtained from Eq. 
solution of ( 40 ) as a; — )■ 00 if the asymptotics at a; - 



Q]. The problem is now to find the asymptotic expansion of the 
-)• are given by ( 42 ) . This problem was addressed in Ref . i44j who 



argued that the large-x asymptotics may be found in the form (31|, (32 1 [i5] . 

It is then straightforward to calculate the coefficients Rn,m of the expansion (31 ), (32) order by order, by substituting 
this expansion into the Painleve V equation ( 40 1 and starting from the first two known terms of the expansion [461 . 
At each order, a second-order differential equation is solved for /„(K,a;), with the integration constants fixed by 
requiring that no extra terms are generated beyond those in Eqs. (|31|), (32). After that, the derivation fully repeats 
that in Section 



calculated coefficients R 



IV we can restore the coefficients /«(«;) in our 
by using relations (37). 



'periodic" form of the expansion illh (12 1 from the 



Of course, this method produces the same asymptotic expansion 
as the method of Section IV We have checked iFoy an explicit calculation of the coefficients up to n = 7 and 

found that the two methods give identical results (14). Furthermore, we extended the calculation via the Painleve V 
equation up to n = 15: this allowed us to calculate the first 15 coeffiencts /„(«) and verify the expansion (II|-(12) 
up to x^^^. 



VI. SUMMARY AND DISCUSSION 



The main result of this work is the conjecture of the "periodic form" of the asymptotic expansion (1I)-(I2| for 
the Wiener-Hopf determinant ([8|. This expansion is based on the proven form p4[ ), together with the "periodicity 
conjecture" (35). The latter can be verified to any order in 1/x by an explicit calculation. 



Our methods only allow to establish the asymptotic expansion away from the Fisher-Hartwig switching points 
= ^ + 1/2. However, we believe that it also holds there [this is supported by the numerical calulcation of Ref. [5] 
in the lattice case] , which would result in a uniform (in k) estimate: 



x{k, x) ^ ^ exp 

j = -imax 



2i{K + j)x ~ 2{k + jf \nx + C{k + j) + £ /„(k + j) {ixY 



71=1 



where 



2 2 



+ o(x-"°), (43) 



(44) 



[•] denotes the integer part, and we assume — f/2 < k < 1/2. 

In the context of the generalized Fisher-Hartwig conjecture, our expansion ( [TT| ), (12) may be viewed as a detailed 
description of the switching between Fisher-Hartwig branches as a function of the parameter n. In physical terms, 
this switching between asymptotic branches is a particular example of a more general notion of "counting phase 
transition" introduced for full-counting-statistics problems in Ref. ITfl 

A practical application of our result is a convenient method of computing cumulants of the number of fermions 
Q on a line segment of length L in the free-fermion problem described in Section n] While those cumulants may, in 
principle, be computed using the Wick theorem (see, e.g., Ref. [S|, the complexity of such calculations grows rapidly 
with the order of the cumulant and with the degree of the 1/L correction. Remarkably, the same cumulants may 
be obtained by the straightforward Taylor expansion of our series (III, (12) at k = 0. Below we list first several 
cumulants up to the order x~^ obtained in such a way: 



-'HQ')) 
-'HQ')) 
-'HQ')) 
-HiQ')) 
-'HQ')) 



X, 

1 + A 
3 3 



— 5" cos(2a;) 



1 



2x 2x2 
3.,„. 15 
2''"' 4x 



2_^3 sin(2x) -I- o{x 3) , 

1 3 , 

Asin(2a;) - ^ + t^(3 - 4A) cos(2x) -I- o{x~^) 



8x3 42;3 



-oC(3)-,^, 

5 r . . 165 5 



'-A^ cos(2x) - ^(3 -t- 4A)(3 - 2A) sin(2x) -t- o(x-^) 
^ 2x-^ 



(45) 



[C(3) + 8A3] sin(2. 



X) + — + — [C(3) - A(9 + 2A)(9 - 4A)] cos(2x) -t- o(x-^) , 
ox-^ x-^ 



2x2 

^C(5) - ^A [C(3) + 2A2] cos(2x) + ^ [(3 - 4A)C(3) - 2A\'S - 2Af] sin(2x) + o(x-3) , 



where 



A = log(2x) + -fE , 



(46) 
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X = kpL as before, is the Euler-Mascheroni constant, and C(n) is the Riemann zeta function [which arise from the 



expansion of the Barnes G function in Eq. ( 13 1] 
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